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ABSTRACT 

Growing evidence indicate supermassive black holes (SMBHs) in a mass range of 
-^BH~ 10^ — lO^^M© lurking in central stellar bulges of galaxies. Extensive observa- 
tions reveal fairly tight power laws of Mbh versus the mean stellar velocity dispersion 
a of the host stellar bulge. Together with evidence for correlations between Mbh and 
other properties of host bulges, the dynamic evolution of a bulge and the formation of 
a central SMBH should be linked. In this Letter, we reproduce the empirical Mbb. — <J 
power laws based on our recent theoretical analyses (Lou & Wang; Wang & Lou; Lou, 
Jiang & Jin) for a self-similar general polytropic quasi-static dynamic evolution of 
bulges with self-gravity and spherical symmetry and present a sensible criterion of 
forming a central SMBH. The key result is Mrh = cr^/ti-") where 2/3 < n < 1 
and £ is a proportional coefficient characteristic of different classes of host bulges. By 
fitting and comparing several empirical A/bh ^ c power laws, we conclude that SMBHs 
and galactic bulges grow and evolve in a coeval manner and most likely there exist sev- 
eral classes of galactic bulge systems in quasi-static self-similar evolution and that to 
mix them together can lead to an unrealistic fitting. Based on our bulge-SMBH model, 
we provide explanations for intrinsic scatter in the relation and a unified scenario for 
the formation and evolution of SMBHs in different classes of host bulges. 
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1 INTRODUCTION 

It is now widely accepted that supermassive black holes 
(SMBHs) within a mass range of 10^ ~ 10^° Mq (M© = 

2 X 10"^^ g is the solar mass) form at the centres of spi- 
ral and elliptical galaxies (e.g., Lynden-Bell 1969; Kor- 
mendy & Richstone 1995; Kormendy 2004). Observation- 
ally, SMBH masses Mbh correlate with various properties 
of spiral galaxy bulges or elliptical galaxies, including bulge 
luminosities (e.g., Kormendy & Richstone 1995; Magorrian 
et al. 1998; Marconi & Hunt 2003), stellar bulge masses 
Mbuige (e.g., Magorrian et al. 1998; Marconi & Hunt 2003; 
Haring & Rix 2004), Mbh versus Mbuige relations in ac- 
tive (AGN) and inactive galaxies (e.g., Wandel 1999, 2002; 
McLure & Dunlop 2002), galaxy light concentrations (e.g., 
Graham et al. 2001), the Sersic index (Sersic 1968) of sur- 
face brightness profile (e.g., Graham & Driver 2007), inner 
core radii (e.g., Lauer et al. 2007), spiral arm pitch angles 
(e.g., Seigar et al. 2008), bulge gravitational binding energies 
(e.g., AUer & Richstone 2007) and mean stellar velocity dis- 
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persions cr (e.g., Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Tremaine et al. 2002; Ferrarese & Ford 2005; Hu 2008). 
These empirical correlations strongly suggest a physical link 
between SMBHs and their host bulges (e.g., Springel, Mat- 
teo & Hernquist 2005; Li, Haiman & MacLow 2007) . 

Among these empirical relations, Mbh and a correlate 
tightly in a power law with an intrinsic scatter of < 0.3 dex 
(e.g., Novak et al. 2006). This relation was explored theoret- 
ically (e.g.. Silk & Rees 1998; Fabian 1999; Blandford 1999) 
before observations (e.g., Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000) and the emphasis was on outflow effects for 
galaxies. The idea was further elaborated by King (2003). A 
model of singular isothermal sphere with rotation (Adams, 
Graff & Richstone 2001) was proposed for the Mbh — o" 
relation. This relation was also studied in a semi-analytic 
model (Kauffmann & Haehnelt 2000) with starbursts while 
SMBHs being formed and fueled during major mergers. Ac- 
cretion of collisional dark matter onto SMBHs may also give 
the Mbh -o- relation (e.g., Ostriker 2000; see Haehnelt 2004 
for a review) . There are also numerical simulations to model 
feedbacks from SMBHs and stars on host galaxies. 

Observationally, there are two empirical types of bulges: 
classical bulges (spiral galaxies with classical bulges or ellip- 
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tical galaxies) and pseudobulges (e.g., Kormendy et al. 2004; 
Drory & Fisher 2007). While SMBHs in classical bulges are 
formed after major mergers, pseudobulges do not show obvi- 
ous merger signatures. Interestingly, pseudobulges also man- 
ifest a Mbh — o" power law yet with a different exponent (e.g., 
Kormendy & Gebhardt 2001; Hu 2008). 

Self-similar dynamics of a conventional polytropic gas 
sphere has been studied earlier (e.g., Yahil 1983; Suto & Silk 
1988; Lou & Wang 2006 (LW06) and references therein). 
LW06 obtained novel self-similar quasi-static dynamic so- 
lutions which approach singular polytropic spheres (SPS) 
after a long time. Such polytropic dynamic solutions have 
been further generalized and applied to study protostar for- 
mation, "champagne flows" in H II regions, stellar core col- 
lapse, rebound shocks and the formation of compact stellar 
objects in a single fluid model (LW06; Lou & Wang 2007; 
Wang & Lou 2007, 2008; Hu & Lou 2008) as well as galaxy 
clusters in a two-fluid model (Lou, Jiang & Jin 2008; LJJ 
hereafter). Here, we take such quasi-static solutions in a gen- 
eral polytropic fluid to describe long-time evolution of host 
stellar bulges and the formation of central SMBHs and to 
establish Mbh — cr power laws. 



2 A POLYTROPIC SELF-SIMILAR DYNAMIC 
MODEL FOR A/bh a POWER LAWS 

For the dynamic evolution of a stellar bulge in a galaxy, 
we adopt a few simplifying assumptions. First, we treat the 
stellar bulge as a spherical polytropic fluid as the typical age 
~ 10^ yr of galactic bulges is long (e.g., Frogel 1998; Gnedin, 
Norman & Ostriker 1999) that they are continuously ad- 
justed and relaxed. Stellar velocity dispersions produce an 
effective pressure P against the self-gravity as in the Jeans 
equation (e.g., Binney & Tremaine 1987). Secondly, the total 
mass of the interstellar medium in a galaxy is ~ lO'^ — 10* M© 
(e.g., Gnedin et al. 1999), only 10"^ ~ 10"^ of the total 
bulge mass. Although gas densities in broad and narrow line 
regions of AGNs are high, the filling factor is usually small 
(~ 10~^; Osterbrock & Ferland 2006) and the gas breaks 
into clumpy clouds. Thus gas is merged into our stellar fluid 
as an approximation. Thirdly, the diameter of broad line re- 
gions of AGNs is only ~ 0.1 pc (Osterbrock & Ferland 2006) 
and the disc around a SMBH is even smaller while a galactic 
bulge size is > 1 kpc. We thus ignore small-scale structures 
around the central SMBH of a spherical bulge. Finally, as 
rotation curves of galaxies show (e.g., Binney & Tremaine 
1987), the effect of dark matter halo in the innermost region 
(around several kpcs ) of a galaxy may be neglected. So the 
dark matter is not included as we study the bulge dynamics. 

Hydrodynamic equations of a general polytropic bulge 
model with spherical symmetry are mass conservation 

dM dM ^ 
_+u—=0 and 



and 'specific entropy' conservation along streamlines (LJJ) 



dM ^ 2 
— — = iTTr p 
or 



(1) 



or equivalently 



dp 1 



dt 



dr 







(2) 



radial momentum conservation (LW06; Wang & Lou 2008) 

1 dP _ GM 
p dr 7-2 ' 



du du 
dt dr 



(3) 



d d \ f P\ „ 



(4) 



where r is radius and t is time; M(r, t) is the enclosed mass 
and u{r, i) is the bulk radial flow speed; P(r, t) is the effective 
pressure and p(r, i) is the mass density; 7 is the polytropic 
index for the stellar bulge fluid; G is the gravity constant. 

As the bulk flow of stellar fluid is slow, we invoke 
the novel self-similar quasi-static solutions (LW06; LJJ) to 
model the bulge evolution under spherical symmetry. We 
use a self-similar transformation (LW06; LJJ) to solve the 
general polytropic fluid equations (O — Q, namely 



r^l/2 ,n 

r = K ' t X 



a[x) 



P = 



AtvG 



M = 



(3n - 2)G 



(5) 



Here, x is the independent dimensionless similarity variable 
while K and n are two scaling indices^ ct[x) is the reduced 
mass density and v{x) is the reduced radial flow speed; P{x) 
is the reduced pressure and m{x) is the reduced enclosed 
mass; reduced variables a, [3, v and m are functions of x 
only. We require n > 2/3 for a positive mass. 

By transformation ([5]), we readily construct self-similar 
quasi-static dynamic solutions from the general polytropic 
fluid equations (HI ~ Q that approach the SPS as the leading 
term for small x. Properties of such asymptotic solutions 
to the leading order (LW06; Wang & Lou 2008; LJJ) are 
summarized below. Both initial {t 0^) and final [t — > 
-l-oo) mass density profiles scale as ~ r~^^"; accordingly, 
the bulge enclosed mass profile is M oc r^~^'^". As r ^ 0^ 
or t ^ -l-co, the density and enclosed mass profiles are 



M = 



(3n - 2)G 



(3n-2)/n 



p — —K ' r 



2/n 



(6) 



where A = {n'^-" /[2{2 - n)(3n - 2)]}i/(«+^-2) and q = 
2(n -I- 7 - 2)/(3n - 2). For either x ^ 0+ or x +00, 
the reduced velocity v ^ 0, which means at a time t, for 
either r — > 0"*" or r — > +00 the flow speed — > 0, or at 
a radius r, when t is short or long enough, the radial flow 
speed u ^ 0. Our model describes a self-similar bulge evolu- 
tion towards a nearly static configuration after a long time 
lapse, appropriate for galactic bulges at the present epoch. 

As the effective pressure P results from stellar bulge 
velocity dispersion, we readily derive the mean velocity dis- 
persion a in a bulge. By specific entropy conservation along 
streamlines, we relate P with p and M (LJJ) and derive the 
P profile from our quasi-static solutions. We then take the 
local stellar velocity dispersion as aL{r,t) = {'yP/ p)^^^ . The 
asymptotic expression when t — > -f 00 for the local stellar 
velocity dispersion in our model is therefore 



CTL 



(j.) ^ ^l/2^1/{2n)^q/2^{9+7-l)/2^{n-l)/„ 



(7) 



To compare with observations, we derive the spatially av- 
eraged stellar velocity dispersion a within the bulge. The 
bulge boundary is taken as the radius rc where p drops to a 
value pc indistinguishable from the environment. Practically, 
the mean velocity dispersion a is usually estimated within 

^ Here, n is not the Sersic index of surface brightness profile. 
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Figure 1. The criterion of forming a central SMBH in a self- 
similar quasi-static bulge evolution in a general polytropic formu- 
lation (LW06; LJJ). The enclosed mass power law is M oc 
with n = 0.751 (solid curve). Meanwhile, we draw a straight 
dashed line A/g = frc^/(2G) with / = 1 for the mass of 
the central SMBH versus the Schwarzschild radius r. Here at 
Tb = 7.2 X 10^ km, the straight line intersects the enclosed mass 
power-law curve for a 2.45 X IO^Mq SMBH. The region within 
rs = 7.2 X lO'^ km is the SMBH in this example. A self-similar 
general polytropic quasi-static solution of n < 1 can thus form a 
SMBH at the stellar bulge centre by this criterion. 



the half-light radius which is less than the outer bulge photo- 
metric radius. As ctl decreases with increasing r by equation 
([7]), the mean a would be higher within the half-light radius. 
In principle, the environmental density for different bulges 
is not the same and it is difficult to give the density exactly. 
Here, we take a reasonable critical density pc as the criterion 
to define the boundary of a bulge, which means when the 
density of a bulge drops to the critical density pc, we define 
the domain within this radius rc as the bulge. For a class of 
bulges with the same n, pc is regarded as a constant. So the 
difference in pc refiects different environments of bulges with 
different parameter n. One can readily show that within rc 



47rr3 



CT_L(r)47rr dr 



)/2^3n9/4^1/2 



(8) 



This relation is particularly satisfying on the intuitive 
ground. A SMBH forms at the centre of a galactic bulge that 
evolves in a self-similar quasi-static manner. Such a SMBH 
was formed by the core collapse of collections of stars and 
gas towards the bulge centre and grows rapidly by matter 
accretion at an earlier phase (e.g., Lynden-Bell 1969; Hu et 
al. 2006; see Haehnelt 2004 and references therein for a re- 
view of the joint formation of SMBHs and galaxies) . As the 
growth timescale for SMBHs is only ~ 10^ yr, our quasi- 
static solutions describe the relatively quiescent phase of 
galactic bulges after the formation of SMBHs as a longer 
history of a bulge evolution. The stellar fluid made up of 
stars and condensed gas clouds has a slow bulk flow speed 
towards the central SMBH, sustaining a reservoir of mass 
accretion for the circumnuclear torus and/or disc on smaller 



scales. Besides, there have been observational evidences to 
show that stars torn up by the tidal force of the SMBH may 
account for the observed X-rays at the centre of bulges (e.g., 
Zhao, Haehnelt & Rees 2002; Komossa 2004; Komossa et al. 
2004; Komossa et al. 2008), which presents a scenario of how 
our stellar fluid being accreted by a central SMBH. 

We now introduce the criterion of forming a SMBH 
in a spherical stellar bulge. A SMBH mass Mbh and its 
Schwarzschild radius are related by Mbh = /r,cV(2G) 
where c is the speed of light and / is an adjustable fac- 
tor of order unity. According to equation ([6]), the mass en- 
closed within the radiufl r is Af oc r^~^^"-. At a certain 
radius r = = [(3n - 2)fc^/(2nAK^^'')]"^^'^"-'^^ where 
it happenflfl M = rscV(2G), a SMBH forms. Only those 
quasi-static bulges with n < 1 can form central SMBHs 
as shown in Figure [T] we derive A/bh = [fc^/{2G)][{3n - 
2)/cV(2nA)]"/(2"-2'A'^/(2-2") and the power law below 



Mbh 



nA 



3n~2 



n/(2-2n) 



2 

7^ 



l/{n-l) 



l/{l-n) 



(3n-2)/(2-2n) 



l/{l-n) 



G 



(9) 



where the coefficient C depends on fc , G, n, 7, pc, and the 
exponent 1/(1 — n) > 3 because of the requirement 2/3 < 
n < 1. Other researchers have their own criteria, which differ 
from ours, to define the black holes in their models to give 
a Mbh — o" relation (e.g., Adams et al. 2001). 

While the Mbh — o" power law is very tight with intrinsic 
scatter ^ 0.3 dex for SMBHs and host bulges (e.g., Novak & 
Faber 2006), such scatter is large enough to accommodate 
different exponents in the Mbh — o" relation (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000; Tremaine et al. 2002). 
By relation ^ , we have a natural interpretation for intrinsic 
scatters in the observed Mbh — o" power law. In our model, 
all bulges with the same n lie on a straight line with the 
exponent 1/(1 — n) as shown in Figure (2] For a fixed n, dif- 
ferent bulges are represented by different K values in trans- 
formation (O, leading to different Mbh and a. However, for 
bulges with different n values, they lie on different lines. For 
elliptical galaxies or bulges in spiral galaxies, they appear 
to eventually take the self-similar evolution described above 
with a certain n value. But pseudobulges may take on differ- 
ent n values. Observationally, we cannot determine a priori 
the specific n value for a bulge but simply attempt to fit all 
bulges with a single exponent, which then leads in part to 
intrinsic scatter in the observed Mbh — o" power law. 

To show this, we fit three published Mbh — o" 
power laws in Figure [2] The first one is Mbh = 
1.2 X lO*M0(a/2OO km s-i)3-^'5 given 
al. (2000) with our parameters {n 
{0.733, 1.327, 0.47 A^opc"^}; and the five points (aster- 
isks *) correspond to K = {0.8, 1, 2, 3, 4} x 10^^ cgs 
unit and rc = 0.73, 0.82, 1.15, 1.41, 1.63 kpc. The sec- 
ond one is log(MBH/M0) - 
given in Tremaine et al. 



in Gebhardt et 
7, Pc} being 



8.13 + 4.021og(cr/200 km s~^) 
(2002) with our parameters 



^ In the domain of bulges and SMBHs, these asymptotic expres- 
sions represent a very good approximation of the exact solution. 
^ We may take the last stable orbit as the cutoff radius or a mildly 
relativistic case of Mbh = frsC^/{2G) with 0.2 < / < 3 to form 
central SMBHs, and this would not alter our results significantly. 
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Figure 2. Power-law Mbh — cr relations in our general poly- 
tropic model for quasi-static self-similar bulge evolution and the 
formation of a central SMBH. Three indices n = 0.7537, 0.7512, 
0.7333 are adopted for three different classes, namely, solid line 
(Hu 2008), dashed line (Tremaine et al. 2002), dotted line (Geb- 
hardt et al. 2000), of Mbh— f power law l[9]l. Given n, we calculate 
the SMBH mass Mbh a-nd the moan stellar velocity dispersion a 
for a certain K value in self-similar transformation JSj. A bulge 
has different mean velocity dispersion a for different K values. 
Each bulge-SMBH system is represented by a point in this plot. 
All such systems of same n lie on a straight line, while systems 
of different n values correspond to different lines in our model. 



{n, 7, pc} being {0.7512, 1.330, 0.0122 A^Qpc"^}; and 
the seven points (plus signs correspond to K = 
{1, 2, 4, 6, 8, 10, 20} X 10^^ cgs unit and Tc = 
2.92, 4.13, 5.83, 7.14, 8.25, 9.22, 13.04 kpc. The third 
one is log(MBH/M0) = 8.28 4.06 log(cr/200 km s"^) 
given in Hu (2008) with our parameters {n, 7, p^} being 
{0.7537, 1.332, 0.00364 Mqpc"^}; and the six points (solid 
dots) correspond to A" = {0.6, 0.9, 1.5, 2.5, 3.5, 4.5} x 10^^ 
cgs unit and Vc = 5.74, 7.03, 9.08, 11.72, 13.87, 15.73 kpc. 
Clearly, to fit all these points in Figure[2]with a single power 
law, we would get a different result with higher intrinsic scat- 
ter. In fact, there is yet another Mbh — cr relation given in 
Ferrarese & Ford (2005), namely log(MBH/M0) = 8.22 + 
4.861og(cr/200 kms-i) (also Ferrarese & Merritt 2000). If 
we were to fit this relation according to our model (e.g., 
n = 0.7942, pc = 1.46 x lO'^Mopc"^, = 3.25 Mpc), the 
critical mass density p^ would be too small and the bulge 
size for a certain average velocity dispersion would be too 
large. So it seems that our model favours multiple power laws 
contained in the available data. In the three more sensible 
fitting examples above, bulge infiow speeds of stellar fluid 
are slow (~ 0.1 — 1 km s~^), an evolution feature of our 
self-similar quasi-static solutions. Near the SMBH bound- 
ary rs, the infiow rest mass-energy flux falls in the range 
of 10*" - 10*^ erg s ^ in these examples, sufficient to sup- 
ply the observed X-ray luminosities (Komossa et al. 2008). 
There can be outgoing accretion shocks around a SMBH in 
these infiows. As the age of galactic bulges is so long (~ lO' 
yr) (Frogel 1998) that such shocks should have already gone 
outside bulges and dispersed or merged into surroundings. 

Besides the Mbh — cr relation, observations reveal 
Mbh cc Af,J^^g^, with A/buigc being the stellar bulge mass 



(e.g., Haring & Rix 2004 and references therein) and Mbh oc 
E°'^ with Eg being the absolute value of the bulge gravita- 
tional binding energy (e.g., AUer & Richstone 2007). Us- 
ing our criterion of forming a SMBH and the bulge ra- 
dius Tc, we derive a power law between Mbh and Mbuigc as 

0.75, 



A/bh <x Mi^/,!^^ according to equation ((6|. For n 



l/(3-3n) 
bulge 

our result leads to relations in Adams et al. (2001) but for a 
nonisothermal general SPS. The bulge gravitational binding 
energy, without contributions from dark matter halo and a 
disc as in AUer & Richstone (2007), is Eg « /J " GMpA-Krdr. 
For self-similar quasi-static dynamic solutions in general 
polytropic fluid model, we obtain Mbh oc Ey^^~^"\ 

As another class of bulges, pseudobulges are thought to 
have formed without merging in contrast to classical bulges. 
Pseudobulges again follow a Mbh — cr power law (Hu 2008) , 
i.e., log(A'fBH/A/0) = 7.5 4.51og(a/200 km s-i). Pseu- 
dobulges may take a different self-similar quasi-static evo- 
lution for a different n. For their different formation his- 
tory, they show a different Mbh — cr power law as observed. 
For {n, 7, pc} being {0.7778, 1.34, 0.000426 Mopc"^}, 
we can fit the empirical power law for pseudobulges; e.g., 
with _K' = 6 X 10^^ cgs unit, we have a = 124 km s~^ and 
= 5.742 kpc for elliptical galaxies. 



3 CONCLUSIONS AND DISCUSSION 

On the basis of a self-similar quasi-static dynamic evolu- 
tion of a general polytropic sphere (LW06, LJJ), we es- 
tablish Mbh = C a i/ti""' with 2/3 < n < 1 by equa- 
tion ([9]). Our first conclusion is 1/(1 — n) > 3 which ap- 
pears consistent with observations so far. Secondly, uncer- 
tainties in the formation criterion of a SMBH [i.e., factor / 
in A/bh ~ frs(?/{2G)\ and in the choice of Tc and thus pc 
will not change the form of equation ^ and the n value but 
will only affect the value of C. Thirdly, the tight Mbh — cr 
power laws and other relations among the SMBH mass Mbh 
and known properties of host stellar bulges strongly suggest 
coeval growths of SMBHs and galactic bulges (e.g.. Page et 
al. 2001; Haehnelt 2004; Kauffmann & Haehnelt 2000; Hu 
et al. 2006). Fourthly in our model, while forming a SMBH 
at the bulge centre (e.g., by core collapse of gas and stars or 
by merging), the spherical general polytropic bulge evolves 
in a self-similar quasi-static phase for a long time. We then 
reproduce well-established empirical A/bh — f power laws. 
Different energetic processes appear to give rise to different 
scaling index n values, which finally determines the slope of 
the A/bh — cr relations in a logarithmic presentation. 

Besides classical bulges and pseudobulges, there are also 
'core' elliptical galaxies (i.e., those with apparent 'cores' 
of relatively flat brightness; Lauer et al. 1995; Hu 2008), 
thought to have formed by 'dry' mergers (i.e., almost with- 
out gas). A steeper Mbh — cr relation exists in these galaxies 
as compared to that for classical bulges (e.g., Lauer et al. 
1995; Laine et al. 2003; Lauer et al. 2007). 

This can also be accommodated in our unifled scenario 
that all hosts of SMBHs may finally evolve into self-similar 
quasi-static phase with different scaling parameters (i.e., dif- 
ferent index n for the slope and different pc for the normal- 
ization of the Mbh — cr relation). While these bulges can 
be of quite different kinds in galaxies, they have these simi- 
lar tight relations and we thus provide a unified self-similar 
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dynamic framework to model the relatively quiescent evolu- 
tion phase of SMBH host bulges and the growth of SMBH 
masses. As the observed Mbh — f relation for classical bulges 
is tight, the elliptical galaxies and spiral galaxies appear to 
take on close n values for merging processes. 

In our model, n is a key scaling index to determine the 
exponent of the Mbh — o" power law. The smaller the value 
of n is, the steeper the density profile is and the smaller 
the index of the Mbh — cr relation is. If SMBHs are formed 
by collapse of stars and gas and a less steeper density dis- 
tribution may provide a more effective mechanism to form 
SMBHs, then we conclude that for a certain value of veloc- 
ity dispersions, the smaller the mass of an initially formed 
SMBH is, the smaller the value of n is. 

After a long lapse, our quasi-static solutions approach 
the static SPS solution as the leading term that is inde- 
pendent of the timescale. So all the described relations here 
are nearly independent of time, as long as the systems have 
evolved for a long enough time. It is not obvious to decide 
when the host bulges began to take the described self-similar 
evolution. But even if we take different times in our model, 
our results would remain largely the same. The only dif- 
ference is that for an early time, we may have a chance to 
observe accretion shocks within the region of stellar bulges. 
Such shocks are characterized by a rapid inner density rise 
of several times and a rapid inner rise of stellar velocity dis- 
persions depending on the strength of accretion shocks. 
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